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We develop a dynamical approach based on the Schwinger-Keldysh formalism to derive a field- 
theoretic description of disordered and interacting electron systems. We calculate within this for- 
malism the perturbative RG equations for interacting electrons expanded around a diffusive Fermi 
liquid fixed point, as obtained originally by Finkelstein using replicas. The major simplifying fea- 
ture of this approach, as compared to Finkelstein's is that instead of N — » replicas, we only 
need to consider N = 2 species. We compare the dynamical Schwinger-Keldysh approach and the 
replica methods, and we present a simple and pedagogical RG procedure to obtain Finkelstein's RG 
equations. 
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I. INTRODUCTION 

Since the original idea of an impurity driven metal- insulator transition (MIT) was put forward by Anderson [0, 
a substantial amount of work has been carried out to understand this problem in the language of phase transitions 
and of the Renormalization Group (RG) [§pB]. Simple scaling arguments were made for the case of non-interacting 
electrons. It was demonstrated that there would always be a metal-to-insulator transition in three dimensions as a 
function of the disorder strength, whereas in one and two dimensions even the weakest amount of disorder would 
make the system an insulator at zero temperature B. The lower critical dimension for the MIT is d — 2 for the 
non-interacting electron problem, and e = d — 2 expansions have been carried out to determine the critical exponents 
characterizing the phase transition || . 

These scaling ideas were extended to include the effects of electron-electron interactions by Finkelstein || , and later 
by Castellani et. al. || , who obtained RG equations for the conductance, as well as the singlet and triplet interaction 
coupling constants, starting with a diffusive Fermi liquid fixed point (for a review, see Ref. [Q). This seminal work 
defined a field theoretical language to study the simultaneous presence of interactions and disorder. Unfortunately, 
this approach suffered from being inconclusive - since the coupled renormalization group equations flow to strong 
coupling, away from the perturbative starting point of a diffusive Fermi liquid state - and technically quite involved. 
A recent discussion of Finkelstein's replica theory can be found in Ref. ||. 

The recent experimental discovery by Kravchenko and co-workers || of a MIT in a two-dimensional Si-MOSFET 
and the subsequent discovery of such a transition in other two-dimensional |l^ |l2j electron gas systems has rekindled 
interest in sharpening our understanding of the combined effects of disorder and interactions. In this paper we will 
rederive the RG equations for all marginal perturbations of the diffusive Fermi liquid fixed point (there are no relevant 
perturbations). We will do so using the dynamical Schwinger-Keldysh approach. Our results agree with those obtained 
using the replica method [^) and disorder-averaged perturbation theory [[|. As we will show, the Schwinger-Keldysh 
and the replica approaches to studying the diffusive fermi liquid exhibit a very similar structure. It is important to 
note that, in the replica solution, there are no subtleties involved the N — > replica limit at the perturbative level, 
which reflects the fact that that there is no replica symmetry breaking in the diffusive fermi liquid state. This is the 
underlying reason why the Schwinger-Keldysh and replica solutions, as we will see, look very much alike, with the 
dynamical thermal indices for time-ordered and anti-time-ordered fields behaving as simple bookkeeping devices, just 
as in the replica solution. On the other hand, glassy systems which exhibit replica symmetry breaking, such as the 
recently proposed Wigner glass phase Jl^,|l4| should be rather different and more interesting. We will explore these 
directions in the future. In the present paper, we will limit ourselves to the study of the diffusive fermi liquid, which 
we present below. 

We present the Schwinger-Keldysh approach for disordered and interacting electronic systems in section |J, where 
we show that instead of N — ► replicas, we only need to consider N = 2 species of fields. Using this dynamic 
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approach, we derive in Section [II a non-linear er-model for interacting diffusion modes in the Schwinger-Keldysh 
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formalism. For simplicity, we restrict attention to the 'unitary' case considered in Finkelstein's original paper ||: 
we assume that the electrons interact through a short-range interaction and that they are in a magnetic field which 
affects their orbital motion, but that the Zeeman coupling vanishes (this situation might be realizable in GaAs- 
AlGaAs systems). By straightforward extension of the methods used here, we can also treat the more interesting, 
but more complicated, 'orthogonal' case in which the magnetic field is turned off. In a similar vein, we avoid the 
straightforward but more involved extension to long-range Coulomb interactions. This would include applications 
of our approach to the quantum Hall plateau transitions. We emphasize the similarities and differences with the 
replicated a- model. In section IV, we find the diffusive saddle point and bring the a- model into a form which is 
convenient for perturbation theory about this saddle point. This form is used to derive the Feynman rules, given 
in section |v|, which are needed for one- loop perturbative calculations. In section VI. we use these Feynman rules 
to discuss the one-loop renormalization of this er-model. One of our main purposes here is pedagogy, so we discuss 
the diagrammatics in detail, paying attention to the symmetry factors and minus signs which prove to be crucial in 
determining the physics of the model. As we show, the appare nt co mple xity - due to the large number of diagrams - 
can be offset by a systematic enumeration. Finally, in sections VII and VIII, we discuss the resulting RG equations, 
a constraint on them following from the Ward identity for charge conservation, and their physics. 



II. THE DYNAMICAL SCHWINGER-KELDYSH APPROACH TO DISORDERED AND INTERACTING 

ELECTRONS 

In the Schwinger-Keldysh - or closed time path - formalism [p5|-p7t , a functional integral (see also the formulation 
by Feynman and Vernon, Ref. [jTsj) is constructed for the time-evolution of the vacuum state from t = — oo to t = oo 
and back to t = — oo (the Keldysh contour, shown in Fig. [I[ is just one possible path for this evolution). Such 
evolution of the vacuum brings it back to the initial state, and therefore the vacuum-to-vacuum overlap (or vacuum 
persistence Z) in the closed path formalism is 1. Consequently, the functional integral is automatically normalized 
to Z = 1 for any realization of the disorder potential, so disorder-averaged correlation functions can be calculated 
directly from the disorder-averaged functional integral. The price which must be paid is the doubling of the number 
of the fields in the theory; the second copy of each field propagates backwards in time. 

-CD ► ^ ! + CQ 

FIG. 1. Keldysh contour, going from — oo to oo and back to — oo. 

The application of dynamical approaches, such as the Schwinger-Keldysh formalism, to disordered systems has been 
previously proposed. At the classical level, Martin, Siggia, and Rose |19[ used a dynamical approach and explored the 
independence of the classical generating functional on the disorder. Sompolinsky |2C| ] used this formalism to study the 
mean-field theory for the spin glass problem. Quantum versions of the idea of using dynamics to enforce Z = 1 have 
been proposed by Schuster and Vieira and by Kree p^j. However, these proposals were never stated in a language 
suitable for calculations in disordered systems. A concrete formulation was put forward by Horbach and Schon p3[ , 
who developed a time-path formulation for disordered non-interacting electrons. Recently, Cugliandolo and Lozano 
proposed a closed time path formalism for quantum spin glasses. In work complementary to ours, Andreev and 
Kamenev [ ^5| use the Schwinger-Keldysh formalism to address the issue of gauge invariance in disordered interacting 
electron systems. They discuss the single-particle density of states - which we do not - but do not derive the full set 
of coupled RG equations - which we do, following Finkelstein. 

Here we apply the closed time path formulation, and cast it in a language appropriate to consider the effect of 
interactions. We give an extended discussion of the method in Ref. (^6|, where we discuss, in addition to systems with 
a natural quantum dynamics, a non-trivial extention of the method to systems with no natural quantum or classical 
dynamics, such as disordered, interacting statistical field theories. Below we give a short summary of the method. 

The Keldysh contour in Fig. |l| is just one example of a possible contour for which we can achieve our desired 
objective, namely the absence of a denominator in the correlation functions (i.e., Z = 1). Here we will discuss a 
more general class of contours originally devised for the study of field theories at finite temperatures p^-|29|] . These 
formulations are carried out in real-time, and their objective was to circumvent pitfalls from the strong assumptions 
about analytical continuation of imaginary-time (such as Matsubara's) formulations. 
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FIG. 2. Real-time contour separated into four parts that factorize into separate contributions: Ci U C2 and C3 U C4. 

At finite temperature, the contour of figure [l] is replaced by a contour running from t = — 00 to £ = 00 and thence 
to t = —00 — i/3. A contour of this type is displayed in figure p[ For such a contour, there is an important factorization 
property (26), shown in Ref. |p9fl , for the contributions from each piece of the contour to the functional integral Z: 

g _ 27C1UC2 ^CsUCi ^-j^ 

Only C\ and C2 are important in obtaining physical correlation functions. |30| As a consequence of the factorization 
property, we can obtain these correlation functions from Z ClUC2 which satisfies the normalization property: 

z c lU c 2 = ! ( 2 ) 

Hence, even at finite-temperature, we can work with a partition function normalized to 1. Denoting the fields on the 
upper (Ci) and the lower (C2) pieces of the countour by 

0i (t) = 4>{t) , 2 (t) =(f>(t-ia), (t = real) (3) 

one defines a matrix propagator 

<T o [0 o (t,i)0j(t',a/)] >=-iA ab (t-t',x-x') (4) 

where T c denotes ordering of fields according their position along the contour of Fig. g. The particular choice 
a = (3/2 makes the form of the propagators especially simple. Let <p denote a complex field, bosonic or fermionic. 
The propagator iA ab can be written in terms of the zero-temperature propagator iAg b as 

iA(u;, k) = u{uj) iAo(ui, k) u^u;) , (5) 

where 

iA (u,k)=( iG ^ k) ..JLJ, (6) 



-iGZ{u,k) 
with iGo(uj, k) the usual time-ordered propagator 

^Go(^ - f, x - x') = (T{cf>(t, x)<fl(jf, x')}) (7) 
and —iGQ(uj,k), consequently, the anti-time-ordered one 

-iG* (t-t',x-x') = ((T[<f>(t,x)^{t',x')]))* = (f[(f>(t',x')cf>^t,x)]) (8) 
The matrix u contains the information about the temperature. For bosonic fields, this matrix is given by 

«M=»B M =( C ™ h h ^; tlZ:)' and crfC = T3 ^- (9) 

For a fermion field, the matrix is altered to account for the fermionic statistics, and we have 

u(w) - u F (w) = ( C ° S M " Sin ^ 1 , where A9^ = % - C=°, and cos 2 C = J • (10) 

v ; v ; I— Asm0 w cos, AO u J ' w w ' w i + e -(«-At)/r ^ ' 

Notice that at zero temperature ub,f = !■ 

All correlation functions in a theory with propagator A a b and interaction £[0, <jf \ are obtained using the following 
functional integral over two fields 4>\.2'- 



Z = I D(j>iD(f)2 exp 



1 

2 ra 



(11) 



where all time integrals are to be done from —00 to 00. This gives rise to a crucial relative sign in how the interaction 
comes in for the forward propagating field 4>\ and the backward propagating one <f>2 ■ 
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III. NON-LINEAR SIGMA MODEL WITHIN THE SCHWINGER-KELDYSH FORMALISM 



We start by deriving a non-linear sigma model in the language of the dynamical Schwingcr-Kcldysh formalism. The 
steps are quite similar to the deriva tion in the replica formulation. 

As shown above, we formulate the problem in terms of a real time path integral 



Z = J D^Di/j e iS ^ 



(12) 



for a fcrmionic field 



comprised of two components ipi labelled by the thermal indices i = 1,2 (in addition to other indices such as spin). 
The index i = 1 corresponds to time-ordered fields, and i = 2 to anti-time-ordered ones. 



A. The free action 



The free part of the action for these fermions can be written as 



S [ip\ip] = I d d x I ipi(x,w) u f (lj) cr 3 



w + ^ + e F + if l s S nw 




U! - 





£ + e F - ir, sgnw 



from which we obtain the free matrix propagator 

= iA(u), k) — up(uj) iA (oj, k) u F (u) 

where 



iA (LJ, k) = 
and iGo is the time-ordered propagator 

iGo(u>, k) 



iG (u,k) 

-iG* {w,k) 



(13) 



(14) 



(15) 



Notice that the anti-time-ordered propagator is the adjoint of iGo- The matrix up(uj) containing the temperature 
dependence was defined above. In particular, uf{u) — lthormai at T = 0. In this case we have 



where in the last line we defined r/5 = ip^cr^. 



W + — + €f ] 1 + if] CT3 sgnw 



(16) 



B. Disorder term 



The static disorder contribution to the action is 



Sv[^,i>] = Jd d xj^V(x)^(x,u ) )a 3 iP(x,u ) ) . 



(17) 



Notice that any interaction term entering the action has the following properties: 

1. It does not mix the two fields tpi 2- There will be mixing after disorder averaging, but not before. 
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2. The contributions from the two thermal components i = 1,2 enter with opposite signs (thus the 173) due to the 
negative sign for the time integration along the anti-time-ordered branch. 

The disorder potential is assumed to be Gaussian distributed according to 

P[V(x)]oze-^I dxV2(x) . (18) 

We then integrate out the disorder, using the fact that the Schwinger-Keldysh formulation cancels the denominator 
problem: 

thus generating a four fermion term 

Srandh^,^] = « | / d d X J ^ ift (x, Lj) CT 3 1p(x, Ul) V> f (x, Ul') <7 3 1p(x, Ul') (20) 

/ ^ / 2T~2T ^( x > u ) i>( x > u ') i>(x>u') ■ 



2 

Let us now introduce the Hubbard-Stratonovich matrix field Q that decouples the four fermions (the matrix Q is 
chosen to be Hcrmitian) . We will choose the decoupling in the diffusive particle-hole channel. As in the case of the 
original diffusive Fermi liquid study by Finkelstein, here we will not consider the cooperon contributions. We have 



J DQ e -£;J ddxtl Q 2 e iSHs[Q,#,1>] (21) 



where 



/pi A J 
ddx J 2T u) {x) ^ J) ■ (22) 

The matrix Q has indices in three separate spaces, i.e., it is assembled as a direct product in energy, thermal, and 
spin spaces: 



hcrmitian. 



Thus, the trace of Q corresponds to 

2 



tr Q2 = fdudtS y y ijt ^ ji 0a 

^ J 2n 2tt .4^ ^ 



i,j=l a,/3=l 



When we write Q^' we mean a matrix whose elements are matrices in thermal (i, j = 1, 2) and spin spaces (a,(3 = 
1,2). 

We can then write an expression for Z as a functional integral over ip, ip and Q: 

Z = JdqJ D{f;D^e~^f ddxtlQ2 e iS °^ e^nsW,*,*] (23) 
where, writing So and Shs explicitly in terms of sums over the energy, thermal and spin indices 

So + S H s = J ddx J ^~2~T w ) | ( w + ^ + £f ) S v + ir > CT 3 sgnuJ Saf3 Sojw ' + i( 3^cy /3 ( x )| ^iA x ^ u ') ■ 

(24) 

Upon integrating out the fermions we obtain an effective theory for the matrix Q field, 

Z = J DQ I ddx tr q2 ef ^ trln [~ <n ~ i {^ +€F ) 1 +^+Q(^)] ) ( 2 5) 
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where we introduced f2 and A as matrices in the energy (g) thermal ® spin spaces: 

= w <W S i] <W (26) 

and 

A Lt" /3 = s § n w <W °"3 <W or A = S 3 ® ^3 ®^ (27) 

Energy Thermal Spin 

with £3 defined as matrix in energy space that is diagonal and whose entries are +1 for positive energies and —1 for 
negative ones. 

The expression Eq. (^5|) is the non-linear sigma model for a non-interacting system. The information necessary to 
dispose of disconnected loops generated by the disorder average is encoded in the matrix structure for and A. This 
structure ensures that we can work with only two indices, i,j — 1,2, instead of N replicas with N — » at the end. 
Notice, however, that very much like in the replicated sigma model, the disorder connects Q components sitting on 
different indices. Next, we will add the contributions due to interactions. 



C. Interaction term 



Let us consider the case of short range interactions, in the singlet and triplet channels: 

Sint = Si + S2 

where 

Si[ip\ip]=TvTi J d d x J T^Tj^T 7j~ ^i,oc{ x M ipiA^u + e) of ipt jj3 (x,o/ + e) ip jt0t (x,u') 

= itTi J d d x J 7^7^ 7^ X! ^i,a(x,uj) ip it/3 (x,u + e) of rpj,p(x, u' + e) ipj, a (x,u') 



(28) 



and 

,<?oP>/jt ih] = -TT-Fn / d d T I 

2ir 2?r 2tt 



S 2 [^\tp] = -nT 2 I d d x I ^7^77^ Y ^i,a( x ^ u ) V>i,a(^,w + e) of ip] t p{x,u' + e) ip jt p(x,u') (29) 

ij,a0 



-7iT 2 ( d d x j '^1^7^ X! V'i.a^jW) ip iiCl {x,u) + t) of ' 1pj t p(x,Lj' + 6) 1pj t p(x,U}') 

ij,af3 



Notice that 



1. Once again, the interaction does not mix the fields ip in different thermal indices i ^ j) this is enforced by the 
03 matrix in thermal space. 

2. The contributions from the two thermal components i = 1, 2 enter with opposite signs due to the negative sign 
for the time integration along the anti-time-ordered branch. 

3. Because all four fermions sit on the same index i, using the definition ipj. a = ipj a of does not change the 
structure of the interaction as the positive or negative signs from the 173 matrices always come in pairs. 

We now introduce two Hubbard-Stratonovich fields X and Y to decouple the four fermion interactions. 

giSit&V] = J dy e tSy[Y] e iV27TTl I ddx I J2n, a0 ^("V") ^,0 (*,">') ( 30 ) 

e iS 2 [^,V] _ f DX e iSx[x] e V27rF2 / ddx I S w , a /3 ^.«( x ^) x ^J e ( x ) 1>fA x > u ') ^-n 



where 

X^ = X*.^( W - to') Sa = Y> - <S) S a p S i: 



G 



Notice that the matrices X 1 ^^ and Y^T 13 depend only on the energy difference w — u>' . We choose X Hermitian 
(X = X'), and Y to be anti-Hermitian (Y = — Y T ); we could alternatively, have absorbed a factor of i into Y and 
made it Hermitian. The action for the matrices X and Y is 

S K [X] = ±Jd d xJ^- X^(e) 4 j X^(-e) (32) 

Sy[Y] = \jd d xj ^ Y\e) Y'(-e) (33) 

The action for the i = 1,2 components of X and Y in Eqs. ( |3^j33] ) come with opposite signs, which is necessary 
to generate the correct sign in the four fermion terms for time-ordered and anti-time-ordered pieces. The fact that 
at least one of the components is not positive definite is not a problem here because we work with a real time path 
integral and have the i factor in front of the actions S x and S y as in Eqs. ( ^0| , ^l|) . 

Including the fields X and Y together with the matrix field Q, and integrating out the fermions, we can write the 
partition function (vacuum persistence) as: 

Z = J DQDXDY e -^f ddx tr Q2 e lS ^ e lS ^ 



D. Expansion around the diffusive Fermi liquid fixed point 

We will now expand the trln aroung the Q saddle of a diffusive Fermi liquid. The saddle point solution for the 
non-interacting problem is Q = iruvo A, where vq is the density of states at the fermi level and A is given in Eq. j2^). 
The matrix Q can be rescaled so that the saddle point can be parametrized as Q — WAU. The matrix Q in the 
non-linear sigma model then satisfy the constraints 

Q = Q\ Q 9 = l, trQ = 0. (35) 

The expansion is achieved by allowing a slow position dependence in the matrix U — ► U(x), and considering small 
couplings Ti.2 and small energies lu. Expanding Eq. (|34|) around the saddle point, we have 

DX DY DQ e tS * [x] e lS * [Y] 

e ~i J d d x {D tr (VQ) 2 +4z tr [(ifi-?)A) Q] +2^2tFT^ tr (XQ) +2i\/27rT7 tr (YQ)} ^gg-j 

where we have also absorbed the rescaling of the saddle into the couplings I\2, as well as D and z. Integrating out 
X and Y, we obtain Z = J DQ e~ s ^\ with an effective action for Q as follows: 

S[Q]=S D [Q] + S e -e[Q] 

where 

s D [Q] = lj ddx i D tr ( V ^) 2 + 4z tr K« n - Q]] ( 3? ) 

S e - e [Q] = J d d x [ittT 1 Q 7l Q - ittT 2 Q l2 Q ] (38) 
and the contractions QjiQ and Q72Q correspond to 

1,0.(3 

QI2Q = £ / ^^^^ Q%£ 4 2nS^ c 2 + c 3 - c 4 ) (40) 

At this stage we would like to compare the expressions for the effective non-linear sigma model Eqs. (^7j38). First, 
notice that the structure is very close to Finkelstein's replicated model; however, here we have only two "replicas" 
corresponding to the two thermal indices. Also, the interaction terms QjiQ and Q"f 2 Q do not couple the two thermal 
spaces much in the same way that the interaction does not couple different replicas. In the Schwinger-Keldysh 
approach the two thermal indices appear with an opposite sign, hence the 173*, in contrast to Sij as is the case of 
replicas. 
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IV. PARAMETRIZATION OF THE SADDLE 



We have used above the parametrization Q = WAU, where A = ^Energy ® C3Thermai ® Ispin- The tensor product 
of the S3 in energy space and 173 in thermal space makes the parametrization of the saddle rather cumbersome for 
carrying out calculations. This can be resolved through the use of a transformation that allows us to parametrize 
around the saddle point 



Aw = sgn uj <W Sn S a t 



or A = £ 3 ® g> ^1^ 

Energy Thermal Spin 



In this case the thermal and spin structures for the saddle become trivial, and we need only to worry about the energy 
structure. This transformati on is achieved through 



where 



T 



ij,af3 



Q -> T^QT 



Sij S a f3 i 



The matrix T can be viewed as a direct product of matrices in energy, thermal and spin spaces: T — U(i)® lthermai ® 
Ispin, where U(i) is a rotation in energy space that dependes on the thermal index. This transformation leaves the 
form of the action invariant, with the exception of A — > A and f2 — > fi, where 



n 



uj 5 tl 



?af3 



It is easy to check that the forms of Q"fiQ and Q72Q remain unchanged upon using the transformation T (the 
interaction does not mix thermal indices, so for i = 1 this is trivial, and for i = 2 one can absorb all changes in sign 
by redefining the integration variables). 

The parameterization of Q around the A saddle becomes simple and familiar, as in the work of Hikami BM: 



Q 



J > 



vi - vw 



v 



to' < 

Vi - vv* 



UJ>0 
UJ < 



(41) 



where the matrix V uu > has row indices with negative energies (uj < 0) and column indices with positive energies 
(uj' > 0). As in several other works (see, for example, the work by McKane and Stone p2| , and by Grilli and Sorella 
|p3| ), we will proceed with the study of the nonlinear sigma model by expanding in powers of the V, V* matrices: 

Q = Ee (n) 



where n is the order in the power series expansion in V, V'. For example, 



Q<°> = 



Q 



(i) _ 



W 

V 



)(2) _ 







Q (3) = , 



)(4) _ 







\{VV*)\ 

The quadratic piece of the expansion of the action Sd [Q] in powers of V, V' is 



d d x 
d d x 



D tr (Vg (1) ) 2 + 4z tr [(iCl - r)A) 
duj' 

J .Wo 2tt ^, <0 2tt L J 



2^ 



ITT f- — ' L 



(42) 
(43) 

(44) 
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and 

(45) 
with 

Q V >HQ { ' = > / -7T- -7T- ~0~ ^ ' Q 0,30,4 CT 3 27r5(Wi - W 2 + W 3 - CJ 4 ) (46) 

* — ' J Z7T Z7T Z7T Z7T 

z, a/5 

VZ'Z V£tf °t ~ ^2 + ^3 - W4) (47) 



cr" 27r5(wi - w 2 + ui 3 - ui 4 ) (48) 
V^f, V»£: ^ 2wS(ui - ui 2 + ui 3 - c 4 ) (49) 

Notice that the contribution to quadratic order coming from Q^jQ^ vanishes due to the trace, as do contributions 
to any order in the expansion of Q^fQ where appears as one of the two Q terms. 



,(2) 



[Q] 



d d x 



(i). 



i,ct/3 



Wl ,LJ4>0 



d(jj\ du)4 



dui 2 dui 3 
~2~k^2jx 



Q W 12Q 



(i) 



Edui\ dui 2 dui 3 dui 4 
~2n~2TT ~2jt~2jt 



U)\ ,OJ4>0 



dui\ dui4 



dui 2 du>3 



V. FEYNMAN RULES 



We now present the Feynman rules which follow from the action Sd + S e - e . Our theory is parametrized most 
succinctly in terms of the V matrix used above. The basic propagator is the (VV) propagator or diffusion propagator 
which follows from Sd- It is represented by a thick line, as in figure |^(a), but it can be understood as a double- line 
comprised of an electron and a hole. It carries momentum q and the two frequencies, £, e + w of the electron and hole. 
The arrow is in the direction of the larger frequency. The propagator also carries two thermal indices and two spin 
indices (one each for the electron and the hole). The action So + S e - e contains terms of all orders in V. However, 
not all of these are needed at the one-loop level, as we now show. 

Let us start by counting which vertices we need to include in the renormalization program. We can determine the 
loop order of any diagram ]33| , |34| ] by considering the number n v of vertices with v legs, and the number of external 
lines E. The number of legs in all vertices must equal the number of external legs plus twice the number of internal 
lines (each internal line belongs to two vertices): 

£ + 2/ = ^!jji„. (50) 

V 

The loop order is the net number of momentum integrals, which is the number of internal lines minus the number of 
contraints from momentum conservation in each vertex, plus an overall momentum conservation constraint 

L = I-^n„ + l. (51) 

V 

Thus, 



±v-i)-\e+1. (52) 



With these relations in hand, it is simple to check that we only need to consider three and four point vertices in 
the one-loop (L = 1) renormalization of the sigma model. As mentioned above, the RG flows can be obtained from 
the two-point function alone, so we set E = 2. Thus, we have 



2 

which has only two solutions: (v = 3, n v = 2) or (v = 4, n v = 1) 



1 v-l) (53) 
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Hence, at the one-loop level, we need only consider the two-point vertex which follows from the expansion of S e - e 
to quadratic order in the V matrix; the triangle vertices which follow from its expansion to third order in V; and 
the square vertices which follow from its expansion to fourth order as well as the expansion of Sd to fourth order. 
The index structure of the former follows from (^) . The index structure of the latter two follow from the fact that 
the thick lines are understood as splitting at the vertices to give the thin lines of the vertices; frequencies, spin, and 
thermal indices from the thick lines then split and follow the thin lines, as we demonstrate with the spin indices in 
figure |3|(c). At the interaction 'dots', they follow the same rules as at the two-point vertices. In figure ||(b), there is a 
matrix structure 8 ai 8^ v if the interaction 'dot' is a Ti and <5 aM 5 7l/ if the interaction 'dot' is a T 2 . This is abbreviated 
as r. The frequencies have an overall <5-function at the 'dot'. The Feynman rules are summarized below. 



-3* 




FIG. 3. The propagator and two-, three-, and four-point vertices. 

Propagator and Diffusion Vertex. At the one-loop level, we need only consider the first two terms in the 
expansion of Sd- They are the propagator and the diffusion vertex. The propagator is given by: 

(a) = (V^Z%) V%£(-q)) = 1 ... 5 u V k 5™5^ 2*5^-^) 2nS(u 2 -w 3 ) (54) 

— iz (witTg — <jJ%cj" I + Dq 2 

The second term in the expansion of Sd is the box, (e). It is given by: 

D z r • i 

(e) = — [2 (qi ■ q 3 + q 2 ■ 94) + (<?i + 93) ■ (<?2 + 94)] + ij^ [^1^3 ~ ^a" + perm.J . (55) 

We have omitted in the expression above the momentum and frequency conserving delta functions. The c^'s are the 
momenta on the four external legs. These legs can have arrows alternating between in and out as one goes around 
the box, with two pointing in and two pointing out. 

Interaction Vertices. The two-point vertex results from the lowest order term in the expansion of S e - e : 

(b) = 2m (-Ti 6 a0 6" v + T 2 (5 Q ^) of 2nS (w x - oj 2 + oj 3 - uj 4 ) (56) 

The next terms in the expansion of the interaction are the triangle terms: 

(c),(d) = T«rr (57) 

where T can be either Ti or T 2 (the suppressed index structure is an extension of the two-point case, as we discussed 
above) and the arrowless propagator can have an arrow pointing in either direction. 

Diagrams (f ) (and a similar one with the arrows reversed) and (g) are the quartic terms resulting from the expansion 
of the interaction. 
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(f) = -m^ (g) = -Htt- (58) 

Decoration of the Interaction Vertices. At the one-loop level, we do not need to consider terms of higher order 
than the above. However, there can still be terms of arbitrarily high order in T because there are quadratic interaction 
terms. These can be obtained by grouping all of the quadratic terms together and inverting them to obtain the full 
propagator: 

(V^'2%) V%£(-q)) = DoW 6 a »6^ 2tt5 ( Wl - o, 4 ) 2ttS («a - 

- 2TTiT 1 D 1 D 2 S zj S tl S jk 8 al3 8^ v 2ir5 ( Wl -uj 2 +u 3 - w 4 ) 

+ 2mT 2 D 2 D Q S tj S ll 6 jk 5 av 5^ 2ttS ( Wl - w 2 + w 3 - w 4 ) (59) 

Note that the interaction terms are distinguished from the propagator by the S lJ on the thermal indices and the overal 
5 energy function. Here, z\ — z — 2Y\ + T 2 = z — 2T S and z 2 = z + T 2 = z — 2T t , where we have introduced the singlet 
and triplet interaction amplitudes, T s = T\ — r 2 /2, T t — —Y 2 j2 by rewriting: 

T^S^ - T 2 8 av 5 f3ti = (tx - ir 2 ^ 6^5^ - \ T ^ aP ' ( 60 ) 

The inversion is more simply done in the singlet-triplet decomposition since these matrices are orthogonal. It can 
then be re-expressed in the form (|59|) where D 0l D\, D 2 are given by: 

D„ = - 



iz (cJiCTg — ui 2 a J 3 J j + Dq 2 



-%Z\ (u>iaf — Lo 2 aif^j + Dq 2 



° 2 = 7 " Ti\ " (61) 

— iz 2 ( U)\cr™ — W2CT3 J + Dq 2 

Hence, by inverting the sum of all of the quadratic terms (rather than just the quadratic term coming from So) 
and amputating the propagators entering and leaving the vertex, we effect the replacement 

(b) = 2m ( -ri5j& S al3 S^ + S^S^) af 2it6 (u>i - uj 2 + uj 3 ~ w 4 ) (62) 

V D a Dq J 

This can, equivalently, be obtained by summing the geometric series of figure ^(a) . It is this peculiar property of this 
theory - namely that there are quadratic interactions - which lead to one-loop RG equations which have terms of all 
orders in T 2 . 

(a) 5> + > • > + — #-=3» — • > + • • • 




Similarly, when the interaction in (c), (d), (f), or (g) is on an internal line, we can sum the infinite series of figure 
|](b) or (c) and, again, make the replacement 

r 2 - r 2 -f (63) 

While the decoration of vertices leads to a replacement of propagators ( p3T) , t he decoration of propagators themselves 
- and, in particular, external legs - is superfluous since the first term in(p9L), which has the desired index structure, 
is not affected by decoration. The decoration of an external line of a vertex amounts to propagator decoration and is 
similarly ignored. 

VI. RENORMALIZATION OF THE a-MODEL 

There are four couplings, D, z, T± i2 , in our action, 

S = ^Jd 2 x Dti (VQ (1) ) 2 + tr [(AizCl-rjA) Q (2) ] + J d d x iTiTi Q (1) 7 iQ (1) - iirT 2 Q (1) 7 2 Q (1) (64) 

All of these couplings are present in the quadratic terms in the action. Hence, the RG flows can be computed from 
the two-point function alone, separating the different contributions according to their spin, frequency, and Schwinger- 
Keldysh matrix structures. 

First, there is the diffusion constant, D - or, equivalently, the resistivity, g = l/((2ir) 2 D). In the non-interacting 
(T-model, this is the only coupling constant. In the interacting case, it is still used as the expansion parameter: we 
assume that g is small, but make no assumptions about r^a. The second coupling constant is z, the relative rescaling 
of time versus space, g and z are renormalized by the diagrams of figure ^[ The other couplings are the interactions 
I\2, which are renormalized by the diagrams of figures ^-^3[ Contributions to the Q self-energy determine the flow 
of D or z if they are proportional to q 2 or u, respectively, provided they have the correct thermal index structure. 
They determine the flow of Ti.2 if they are independent of the external and frequency and are diagonal in the thermal 
indices; they contribute to one or the other depending on their spin structure. 

It is a matter of taste whether or not we choose to introduce a 'wavefunction' renormalizatiorfj], £, since this can 
be completely absorbed in the renormalization of D and z. There is a natural definition of £ which is related to the 
density of states, but this is moot for the scaling equations which we consider because we can rescale our couplings 
to eliminate £. Hence, we sweep the wavefunction renormalization under the rug in this paper. 

At tree-level, D, z, T% >2 , are marginal. (One might have expected Tx >2 to have the dimensions of (length) -2 , but 
they do not as a result of the extra frequency integral.) To find the fluctuation corrections to their scaling behavior, 
we use a Wilsonian shell elimination scheme in which we compute their change as we integrate out shells 

A - dA < Dq 2 < A, < lo < A 
and 

A - dA < lj < A, < Dq 2 < A (65) 



1 the quotation marks refer to the fact that we mean a renormalization of V, not ip 
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FIG. 5. The self-energy diagrams which renormalize z and D. 

Renormalization of D and z. Let us examine figure [|. By inspection, we see that these diagrams have the 
appropriate matrix structure for the propagator, i.e. the matrix structure of (p4|). Diagrams without an interaction 
vertex do not contribute to the renormalization of D to one-loop order; the reason is that there is a cancellation between 
the two (time-ordered and anti-time-ordered) Schwinger-Keldysh species. In the case of the replica calculation, the 
contribution to one-loop order is proportional to the replica number N — > 0. 

The diagrams in figure || give the following contribution to the Q-field self-energy, E: 



£(<?, SI + e, e) = E A (g, Si + e, e) + E B (g, SI + e, e) 



where 



« « { -2m) [T 1 (k,u)-2T 2 (k,oj)] D 2 (k,u;) 
(2iry 2n 

is the contribution of the diagram on the left and 
E B (g,Q + e,e)=4 



j(k + q) 2 -i~(u + Sl) 



r d 2 k 


dtu 


f iri\ 


1 (27T) 2 


2^ 





(66) 



(67) 



(68) 



is the contribution of the diagram on the right. There is a relative factor of 2 between the Ti and I^ contributions 
because the former contribution has no free internal spin indices while the latter has one which is summed over. 
Notice that E^ and E# cancel at q = SI = 0. The contribution of E^ as well as a piece of E# are absorbed into 
the wavefunction renormalization. ^] Differentiating with respect to uj or q 2 , respectively, we get, to lowest-order in 
I\2, the one-loop contributions to dD and dz (the external frequencies and momenta have been set to zero after 
differentiation) : 



dD = - 

dz = —i 



/dS 2 



V dq y q=o n= 
fdZ B \ 

J g =o,n=o 



V dSl 



= -D(-2ni) J AL±[r 1 (k,cv)-2r 2 (k ) aj)} Dl{k^)Dk 2 

, L r 1 (fc,o)-2r 2 (fc,o)] D (k,o) 



d 2 k 



(27T) 



(70) 
(71) 



In the dz term, the SI derivative cancels the u> integral. 
Making the replacement (p3j), we get: 



2 The E b contribution includes a term of the form: 

4F2 / ^^ D0{k,Uj)D2{k,LU) m 

which comes from differentiating Do with respect to q 2 . The wavefunction renormalization multiplies the V fields and therefore 
resurfaces in the renormalization of IV 
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dD_ 
~D 



(-2m) J J7^^ Dk2 [r 1 D (fc,w)Di(fc,w)Da(fc,w) - 2r 2 r^(fc, ^(fc, cj)] 

<fa = ~(r!-2r 2 ) I ^ a,(M) 



(72) 
(73) 



The contribution to dz is unchanged but the contribution to dD receives corrections of all order in I\2- This result 
is the full one-loop result: it is lowest order in g but to all orders in 1^ 2- 

The integrals (^Oh and (|72|) are discussed in the appendix. Here, we just state the result: 



dl 

dz 



= g 2 QTi/a (z, Zi, z 2 ) - 2r 2 / 2 (z, z, z 2 )) 



g(-T 1+ 2T 2 ) 



(74) 



where dl = 
functions: 



-d(lnA) and we have made the substitution g 



l 

8jtD 



Following Finkelstein, we have introduced the 



fi (z, z 2 ) = 
h(z,z 1 ,z 2 ) = 



z 2 — z \ z / 



2zi 



Zi - z 2 



-fi (z,zi) 



2z 2 
z\ - z 2 ■ 



-fi (z,z 2 ) 



(75) 



As we show in the appendix, the functions fi. 2 arise as a result of the mismatch between the poles of the propagators 
■Co, i,2- In the non-interacting case, in which these propagators are equal, these functions are just constants, /i (z, z) — 
1, f 2 (z, z, z) = 2. Interactions split the energies of singlet and triplet particle-hole states from their non-interacting 
value, leading to (|75|). 

Renormalization of T\ :2 . The diagrams which renormalize Ti are in figures ^| and |^-p~0| while the related diagrams 
which renormalize T 2 are in figures | and In§|l|, we have, for the sake of completeness and pedagogy, enumerated 
all of the diagrams which contribute to the renormalization of IY2. On combinatoric grounds, we see that we have 
included all 20 of the diagrams which contribute to Ti and the 20 which contribute to T 2 . All of these diagrams 
must be taken into account to obtain the correct signs and relative coefficents of various terms in the RG equations. 
However, these two sets of 20 are really elaborations of two sets of 5 diagrams. In figure we give the 5 basic diagrams 
for Ti; figures |^-Oare put in an appendix to avoid unnecessary clutter and confusion. From diagram |^(a) we obtain 
the 6 diagrams oT||(a). From each of 0(b), 0(c), and 0(d), we obtain 4 diagrams shown, respectively, in figure ||(b), 
|(c), and|(d). Finally, figure 0(e) gives rise to the two diagrams of figure |10|. 

In figure [| Ti contributes to the renormalization of T 2 and vice versa. When only T 2 's appear in diagrams pHl3|, 
the spin index structures dictate that |^[To| contribute to the flow of Ti whi le |lT] - |T^ contribute to the flow of T 2 . If 
any of the vertices are IVs, however, corresponding diagrams in p[]lO| and [Tl|-[13| cancel. Hence, we need only consider 
r 2 's in diagrams 1^11^. 




FIG. 6. The lowest order diagram which renormalizes Ti,2. 
By following the frequencies and spin indices around the diagram of figure ||, we see that it gives a contribution: 

dTiS^S^ - dT 2 5 av 5^ = (T^S^ - T 2 6 al3 6^) f -^4* D Q (k,Lj) (76) 

J (2tt)^ 
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«3- 






FIG. 7. The five 'primitive' diagrams which renormalize Fi. They are elaborated in figures |8|-|10|. The corresponding 
elaboration for F2 is in figures |ll"r[l3| . 



Corresponding diagrams of figures p[]lO] an d]ll| - |l3| give identical contributions, apart from a factor of 2 coming from 
a spin sum, so we will focus on]^-^o|. Figure g(a) gives a contribution 



(77) 



The signs of the six diagrams ((iii) and (iv) are two diagrams each because the arrow can point in either direction on 
the arrowless lines) are determined by the relative signs of the two triangle vertices (p^). The first two diagrams come 
with a positive sign while the other four come with a negative sign. Adding these together, we obtain the coefficient 
—2. Both of the interaction 'dots' are on internal lines, so we make the replacement (p3|): 



2 m dY^^-2^T 2 fJ^^-D\ 



(78) 



In figure ||(b), there are four diagrams and they all come with positive signs. Only one interaction 'dot' is on an 
internal line, so we have: 



(79) 



The diagrams of figure p"l|(b) vanish because the integrated frequency is overdetermined. However, the wavefunction 
renormalization coming from the neglected part of makes up the missing contribution. 

In figure |](c), the four diagrams have positive signs. There are three interactions, but one is on an external line, so 
there's a contribution: 



-.(•O _ a ^2 - I dk 



2m dT^' = 4 (iTTT 2 y (i2nT 2 ) 



A du>i f° duj 2 sl2 
2 / / "TTT D o{k,uJi - uj 2 )[D 2 (k,uJi - uj 2 )\ 



(2tt 



n 27T J- A 2"7T 
A 



(80) 



where we have set lj — loi — lu 2 as the only variable inside the propagators, and obtained the phase space factor ^- 
from integrating the intermediate frequencies. 

Figure 0(d), on the other hand, comes with a negative sign and only one internal interaction: 



2mdT^ = -4 ^T 2 f (z2.r 2 ) / / 



duj LO 
2tt 2vr 



[Da{k,u)] 2 D 2 (k,u) 



(81) 



The coefficient 4 is important for the physics of this theory. In Finkelstein's original paper g, the contributions from 
the diagrams of figure ^ were too small by a factor of 2. As a result, a metallic fixed point with infinite conductivity 
was found. This error was corrected in J36||57]]. 
Finally, we have the two figures |To|(e) 
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2m dT 



« = 2 ^T 2 f (i2nT 2 fj ** g (^) 2 [D Q {k^)f [D 2 (k^)f , (82) 



where, as in cases (c) and (d), there are intermediate frequency integrals that lead to the phase space factor (j^) 2 - 
The integrals in the expressions above are discussed in the appendix. Their upshot is: 

dT< o) =-H s d(lnA) 

Z 2 

dr{ b) = 2T2/i (2, z 2 ) ffd(lnA) 
rflf 5 =-^/ 2 (z,z,z 2 ) gd(lnA) 

z 



dr (e) = _ r , , ^ , d (ln A) (83) 



Adding these contributions, we obtain: 



dr[ a) + dTf ] + dT[ c) + dT[ d) + dr[ e) = -$ = ^ 3 d(lnA) (84) 



VII. RENORMALIZATION GROUP EQUATIONS 



Gathering these results, we have the flow equations 



% = 9 2 ( r i/2 (z, zx,z 2 ) - 2r 2 / 2 (z, z, z 2 )) 



dl 

dz 



= 5 (-r 1 + 2r 2 ) 



di H r2 + f 

dT 2 ( 2Y\ \ 

Adding the second and fourth equations, and subtracting twice the third equation, we find: 

~(«-2ri + r 2 ) = o (86) 

di 

Although we have only verified (^6|) to one-loop, it is, in fact, an exact relation which follows from the Ward identity 
for charge conservation, as we discuss in appendix C. 

Since a constant of integration in ( |S6"| ) can be absorbed into a rescaling of frequencies, we can set Ti = (z + T 2 )/2. 
Changing variables from T 2 to 7 2 = T 2 /z we have the equations: 

d 9 _ 2 (i , o ft 1 + 72 



i=^l + 3^1-^ln( 1 + 72 ; 

dz 



M= Zg \—2 + 2 12 

^=9(l+72) 2 (87) 

^From the final equation, we see that 7 2 increases at long length scales. This is not problematic since we made no 
assumptions about the smallness of T 2 . When 7 2 becomes sufficiently large, z begins to grow, and g, after an initial 
increase, begins to decrease. Since these equations are valid in the small g regime, it would appear that the flow 
improves their validity. However, 7 2 and z diverge at a finite length scale. Since these equations include all orders in 
72, their breakdown signals the onset of non-perturbative physics. 
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It is useful, as a comparison, to consider the BCS interaction in a Fermi liquid. There, the flow equation for the 
BCS interaction, V, can be computed to all orders in perturbation theory: 




The one-loop RG result is the full story because there is only a geometric series of bubble diagrams. This equation 
also breaks down at a finite length scale - the coherence length. At this length scale, the non-perturbative physics 
of pairing takes over. Similarly, the breakdown of our flow equations ( |S5| ) implies that non-perturbative physics - 
such as the formation of local moments - determines the behavior of the system. Such physics cannot be accessed 
perturbatively from the diffusive saddle point. 

VIII. DISCUSSION 

In this paper we have shown how the Schwinger-Keldysh dynamical formulation p(| can be applied in treating 
a disordered interacting electronic problem. We reproduce Finkelstein's Renormalization Group equations for the 
interaction strength and conductance. In our approach, the RG procedure is carried out in a very simple way, in close 
resemblance to a 4 theory (with cubic terms as well). 

For the case of the diffusive Fermi liquid state, the calculation using Schwinger-Keldysh is very similar to the 
replica method. The two thermal indices for time-ordered and anti-time-ordered propagation play a similar role to 
that of the replica indices in the particular case of the diffusive Fermi liquid regime studied by Finkelstein. The 
disorder couples fields with different thermal indices (or replicas), while the interaction does not mix thermal indices 
(or replicas). The Schwinger-Keldysh indices, like the replicas, act in this example as simple book-keeping devices. 
In the replica calculation, the diffusive Fermi liquid state is replica symmetric. This is the underlying physical reason 
for the simplicity of the replica structure. 

There should exist much more interesting saddle points where the replica symmetry is broken. One such point 
is the Wigner glass phase, which occurs in the limit of low electronic densities. As in the cases with other types 
of glasses, such as spin glasses, we know that replica symmetry breaking is connected to broken ergodicity. In the 
Schwinger-Keldysh approach we expect that broken ergodicity (and thus replica symmmetry breaking) will manifest 
itself in subtleties related to the ui — > limit and the symmetry between time-ordered and anti-time-ordered products. 
The Schwinger-Keldysh provides a natural tool to study dynamical effects, which is a direction that we will explore. 
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APPENDIX A: ELABORATION OF THE FIVE BASIC DIAGRAMS 

The five basic diagrams of figure [7] are a schematic representation of the following 20 diagrams which renormalize 

ri. 



17 




18 




FIG. 10. Diagrams (e) renormalizing Ti. 
The corresponding diagrams which renormalize ^ are: 




FIG. 11. Diagrams (a) and (b) renormalizing F2. 
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(iii) >» 





(iii) 





55* •- 



(iv). 




(iv) 





FIG. 12. Diagrams (c) and (d) renormalizing T2- 



(e) 



(i) 




(ii) 




FIG. 13. Diagrams (e) renormalizing F2. 



APPENDIX B: INTEGRALS 



The following five logarithmically divergent integrals arise in the calculation of the diagram of figure |g and diagrams 
(a)-(e). We define x = Dq 2 and integrate over the shells A — dA < x < A,0 < uj < A and A - dA < uj < A, < x < A. 
In the diagram of figure 4, there are no free frequencies; there is only a momentum integration. 



1 



4ttL> 



dx Dq(x,uS) = 



1 



4ttD 



A—dA 



dx — = — !— d(lnA) 



(Bl) 



In diagrams (a)-(e), there are multiple propagators with different z's. By rewriting these products of propagators as 
a sum over poles in x and to, we reduce them to c?(ln A) times non-divergent integrals. These integrals give non-trivial 
functions of the z,-'s. 



1 



2(2tt) 2 D 



dxduj Dq(x, u>) D%{x, ui) 



1 



1 



1 



, dx / duj 

2{2tt)-D JA-dA Jo x ~ i- ZUJ x ~ iz2U> 

1 f A J f A J 1 1 

+ —, — / dbJ I dx 

2(27r)^D J A _ dA J x - izuj x - iz 2 w 
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i , , r A . ( i i 

dA / du 



2(2ir) 2 D J (z 2 — z)lu \A — izuj A — iz 2 u> 



+ 



1 i dA f A , ( 1 1 



/* f— 

Jo \x-i 



2(2ir) 2 D (z 2 -z) A 7 \x - izA x - iz 2 A 

= Inf— ) rf(lnA) 

z 2 - z V z / 2(2tt) 2 L> v ; 

= i/l(z ' Z2) 2(2^D rf(lnA) (B2) 



Similarly, 

If 1 f A f A 1 1 I 

/ x 2n / dx dw xD (x,uj)D 2 {x,uj)D 1 (x,uj) = / dx dw x : : : 

2(27r)^_D J 2(2tt) 2 D J A _ dA J Q x — izbj x — iz 2 bj x - izioj 

r A r A 



1 f , f , 1 1 1 

H ; — — / dtu I dxx 

2(27r)^D J \—dA Jo x ~ i- ZUJ x ~ iz 2 u x — iziLU 



2iz\ , / z 

In 



(-21 - Z 2 ) {Z - Zi) \Zi 

2iz 2 , / z \ \ I 



if2(z,z 1 ,z 2 ) 2 ^l )2D rf(lnA) (B3) 



2{2ttYD J J 2(2t:) 2 D J A _ dA J x - izuj x - iz 2 lu x - iz 2 uj 

r A pA 



1 1 



1 / j r\j 1 1 1 

H t — — / ujdto I dx 

2(2iT) z D J A— dA Jo x ~ i ZUJ x ~ iz 2 oj x — i 



IZ 2 UJ 



2z 2 \ y (z-z 2 ) (z-z 2 ) 2 \z 2 J J 2(2n) 2 D 



2(2 



\ 9 - \ dx dw J 1 D (x, w) D (x, u) D 2 (x, uj) D 2 (x, uj) 



,2 



2(2-k) 2 D J a— dA Jo x ~ x — x — »z 2 cj a; — zz 2 w 

1 /" A 2j / >A J 1 1 1 1 

+ „,„ / oj du) ax 



2(2ir) 2 D JA-dA Jo x ~ x ~ ^ — iz 2 uJ x — iz 2 oj 

(l + l--2h{z,z 2 )\ 1 

= -'( (W J W73 rf(lnA) (B5) 



APPENDIX C: WARD IDENTITY 



In section VI, we found that our one- loop RG equations imply that: 

5= = 4 r - < C1 > 

We would now like to indicate why this constraint follows from charge conservation. Let us imagine returning to 
our disorder-averaged action So + S ran d + ^mt for interacting electrons. Let us consider some fixed realization of the 
disorder. The Ward identity which follows from charge conservation, 
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(T (j„ (y) ^ {xx)i> (a*))) =i8{y~ an) (T (an) V (*a))) - t* (y - a*) (T (tf (an) $ (a*))) (C2) 

relates vertex renormalization (the left-hand-side) to the renormalization of the propagator (the right hand side). 
This is more transparent in momentum space: 



A ; , (;,.<!) = A (a.,- (/1 -v (>)1 ,(-•:-;, 



A M (p, 0) is the F ouri er transform of the correlation function with a current insertion (at zero momentum) on the 
left-hand-side of ( p2| ) . E is the self-energy. We focus on the // = component of this equation: 

Ao(p,0) = ^ («-S(p)) (C4) 

The equality between the left- and right-hand-sides is exemplified by the relation between diagrams |l4|(a) (the dashed 
line indicates the current insertion in the correlation function of the left-hand-side of (p2Q) and fl4|(b). 



As we indicate in figure |l4|(a), the diagrams contributing to the renormalization of the left-hand-side of (C4) are 
equivalent to those which determine the RG flow of T s . 



Meanwhile, the diagrams contributing to the renormalization of the right-hand-side of (C4) renormalize z. One 
way of seeing this is to note that the renormalization of the propagator determines the RG flow of z since the former 
is the renormalization of the term iwi/j'tf) in So while the latter is the renormalization of the equivalent term in Sd, 
namely tr[ujQ]. The factor of 2 on the right-hand-side results from the summation over the spin index a in (b). This 
index is held fixed as an external index in the T s renormalization 

In principle, we can derive the same result for disorder-averaged correlation functions directly within the er-model 
by deriving the corresponding Ward identity. However, this is more cumbersome because the current operator and 
the gauge transformation rules are more complicated in the Q-field language. 
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FIG. 14. Diagrams contributing to the renormalization of (a) T 3 and (b) z. 
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